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Abstract 

We study two correlated electrons in a nearest neighbour tight- binding 
chain, with both on site and nearest neighbour interaction. Both the cases 
of parallel and antiparallel spins are considered. In addition to the free 
electron band for two electrons, there are correlated bands with positive 
or negative energy, depending on wheather the interaction parameters are 
repulsive or attractive. Electrons form bound states, with amplitudes that 
decay exponentially with separation. Conditions for such states to be filled 
at low temperatures are discussed. 
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Exact solutions to problems involving correlated motion of interacting 
particles are extremely rare. Even simple systems like electrons moving in a 
wire are usually solved approximately only. A body of recent literature exists 
for the case of just two electrons moving in a one dimensional disordered 
potential. [1-9] The problem of N particles in an ordered string of length L is 
conceptually simple if one ignores spin, since it is then formally equivalent 
to that of a single particle in N-dimensional space, with the pair interaction 
acting as a defect potential associated to the planes Xi = Xj, where Xi is the 
position of the i-th particle. One then expects a band of about L N extended 
states with finite amplitude in all of space, and bands of L^ N ~ S ^ surface states 
localized about the geometrical defect where S planes intersect, or the planes 
themselves if S — 1. In spite of this qualitative understanding of the ordered 
case, an exact solution has been reported for an N=2 singlet state only. [10]. 
Thus, in this reference, as well as previous work[l-9] it is found that two 
correlated electrons exhibit a density of states characteristic of one particle 
in a plane, that then has a Van Hove singularity at E — 0. 

In this work we report on another exact solution of the two-particle prob- 
lem in an ordered lattice. It describes a paired triplet state, with an energy 
that may fall within the conduction band, making it of interest to supercon- 
ductivity. [11] In order to see how this comes about, we consider a chain of L 
sites within the tight binding model, with up to nearest neighbor interaction. 
With the understanding that the amplitudes c(l,m) represent two electrons 
at sites I and to with either parallel or antiparallel spins, the equation of 
motion without disorder reads, 

—tc(l — 1, 77l) — tc(l + 1, TO) — tc(l, TO + 1) — tc(l, TO — 1) + 

J[Sl+l,m + 8l-l,m]c{l, TO) + U 6i tTn c(l, to) = Ec(l, to) (1) 

A transformation to center of mass coordinates is effected by taking 

c(l, to) = exp(ik(l + m)a)x(l — to). (2) 

Here k is the center of mass momentum and a is the lattice constant. 
Denoting by p = I — to the distance between the two electrons, the equation 
obtained by substituting (2) into (1) is that of a single particle in an effective 
linear chain with sites p, 

-2t cos(ka)x(p - 1) - 2t cos(ka)x(p+l) + J[S p ,i + 6 p -i]x(p) + U 5 pfi x(p) = 
E X (p) (3) 

where U is a contact Hubbard parameter[12] and J the nearest neighbor 
coupling strength. A peculiar feature of this equation is that the effective 
hopping parameter depends on k, actually vanishing at ka = ±|. In the 
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absence of interactions one may set x(p) — exp(ipqa), obtaining the energy- 
band 

E(k, q) = — At cos(ka) cos(qa) (A) 

This covers the range — At < E < At. Assuming L large, Eq.(3) represents 
an a lattice with impurities around the origin and we expect the band (4) to 
hold true even in the interacting case, save for corrections in the density of 
states of order 1/L. One can easily show the dispersion (4) to be exact in the 
case J = 0.[13] 

We will first work out the magnetic case of two parallel spins. The spatial 
wave function must then be antisymmetric under exchange of particles, or 
x(p) = ~x(~P}- We consider equation (3) for the separate cases p — 0, p — 1 
and p > 2, taking x(0) = 0. The equation for p = is trivially satisfied 
because of the antisymmetry of the wave function. In looking for a solution 
for p > 1 we assume there is some constant f3 such that x(p + 1) = Px(p)- 
From the boundary condition that the amplitudes must remain finite as p 
becomes large, one must have | (5 |< 1. Solving the equations we get 

E = J+(^)cos 2 (ka) (5) 

= _(2f ) cos(A;a) (6) 

The solution decays exponentially with exponent 7 = In | J/(2t cos(ka)) \ 
, having the form x{p) = ^pexp(— 7 | p |) where A p vanishes at the origin, 
and has the value A if (3 > and (— if f3 < 0, for finite p. Since 
the amplitude is largest when the particles are next to each other the state 
represents a pair bound state moving with center of mass momentum k. 
Notice that in equations (5) and (6), U is not involved at all, a feature 
already contained in the Hubbard model[14]. Therefore it was necessary to 
include nearest neighbor interaction to find it. Some properties of this paired 
triplet state are the following. First, under a change of sign of the nearest 
neighbor coupling constant, the energy just changes sign. Referring now to 
positive J, at a fixed value of k its energy is above those of the free electron 
band at the same k. The state may still be within such band, yet with 
different center of mass momentum. The lowest energy is J and it occurs at 
I ka I = 7r/2. In this limit state the two electrons are as close as possible, 
with finite amplitude as nearest neighbors only. For other values of k the 
energy is higher, and the pair is larger in size. 

Band overlap between paired and free states occurs for At > J. The 
energy range of overlap is J < E < 2J if 2t > J and J < E < At otherwise. 
The center of mass momentum is bounded by the condition | cos(ka) |< ^ in 
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the former case, and | cos(A;a) |< y |(1 — j^) in the latter. If the Fermi energy 
lies above J, when band overlap occurs, it will be energetically favorable to 
create bound pairs, with wavector around | ka \— | , where the density of 
states has a divergence. The model thus predicts an instability in the Fermi 
liquid with the formation of pairs with parallel spin that coexist with other 
free particle states. Since the pair and the single particle momenta are not 
the same at a given energy, a transfer from one state to the other requires 
some excitation (a phonon for example) to supply the missing momentum. 

Let us go on to consider the case of antiparallel spins. In this case the U 
term is involved, and if (U/J) is large, we can ignore the nearest neighbor 
coupling and solve the model with J = 0, with xip) — x(~P)- Again we try 
x(p + 1) = Px(p) with | (5 |< 1. Following a similar procedure as before we 
find now 

E = ^jU 2 + 16t 2 cos 2 (A;a) (7) 

P = JM (8) 

giving an exponential decay rate 7 = arcsinh|f//4tcos(/ca)|. The sign 
of the square root in equation (7) should be the same as the sign of U . 
This solution was already reported in Ref. 10. Its qualitative features are 
similar to those of the previous one. Again, for positive (negative) U and 
same wavenumber k the values given by equation (7) are above (below) those 
given by equation (4). However, for At > U > 0, allowing for different values 
of k, there may be energies in the paired band lower in energy than some 
in the single electron band. For positive U the lowest energy state for the 
paired band is again at ka = ir/2 and the two electrons are stuck together, 
while as the energy increases above the band minimum U the extent of the 
pair increases as well. 

The most likely case in a real system is U >\ J |> 0. The paired states 
with lowest energy will therefore correspond to parallel spins. Different band 
widths for the paired states are predicted from equations (5) and (7). Effects 
of overscreening[15] may possibly give rise to a positive U, though a negative 
J. As remarked before in the magnetic state U does not enter, and the energies 
for the paired solutions lie below the energies of the single particle band center 
and, for | J \ sufficiently large, even below the conduction band altogether. 

In conclusion, the effect of correlations on interacting electrons, taken two 
at a time, is to form paired states grouped in bands. In contrast to various 
approaches based on Greens functions [17- 18], our treatment is exact, though 
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with the limitation that we only consider correlations between two electrons 
at a time. Correlations among more than two electrons may be important. 
However, based on the equivalent picture of a single electron moving in higher 
dimensions in a lattice with defects suggests that the number of such states 
scales as L~ s . The possibility of having unconventional metals has often been 
traced to the failure of Fermi liquid theory in correlated lower dimensional 
systems. [19-22] It is obvious here that the correlated bands cannot be placed 
in a one to one correspondence with one electron states. Thus these effects 
may also be traced to the possibility of two or more- particle correlations. 
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